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Summary 

The Icing Research Tunnel (IRT) at the NASA Lewis 
Research Center is a subsonic closed-return atmospheric 
tunnel. The tunnel includes a heat exchanger and a refrigeration 
plant to achieve the desired air temperature and a spray system 
to generate the type of icing conditions that would be 
encountered by aircraft. At the present time, the tunnel air 
temperature is controlled by manual adjustment of freon 
refrigerant flow control valves. An upgrade of this facility calls 
for these control valves to be adjusted by an automatic con- 
troller. This report discusses the digital computer simulation 
of the IRT refrigeration plant and the automatic controller that 
was used in the simulation. 


Introduction 

The Icing Research Tunnel (IRT) at the NASA Lewis 
Research Center is a closed-return, atmospheric tunnel 
equipped for testing low-speed aerodynamic models. A stan- 
dard centrifugal compressor refrigeration cycle is used to 
reduce the tunnel air temperature to the desired test conditions. 
At present, the temperature is controlled manually by a 
refrigeration system operator who manipulates the freon 
refrigerant flow control valves. As part of an upgrading and 
rehabilitation of the IRT Facility, this study was performed 
to design an automatic control which could be used as an 
alternative to the manual control. Automatic control would 
make more efficient use of the tunnel run time. Presently it 
takes the system operator considerable time to change the 
tunnel air temperature from one operating point to another. 
This study demonstrates that an automatic controller could 
accomplish this in less time. For this investigation, the IRT 
refrigeration system was simulated on a digital computer. The 
purpose of the simulation was two fold: (1) to test various 
temperature control schemes and design a suitable controller 
and (2) to study future facility automation in areas such as 
startup, shutdown, and decision-making capabilities (such as 
changing the number of compressors during run time). 


System Description 

The major components of the Icing Research Tunnel (IRT) 
Refrigeration System are compressors, pumps, condensers, 
economizers, a flash cooler, and a heat exchanger. The system 
removes heat from the IRT by means of a heat exchanger and 
maintains the tunnel air temperature at the required operating 
point. A standard centrifugal compressor refrigeration cycle 
is used, with Freon- 12 as the refrigerant. The system is operated 
in a flooded manner where subcooled refrigerant is continually 
circulated through the tunnel heat exchanger. This feature 
improves heat transfer by keeping the refrigerant a liquid until 
near the end of its passage through the heat exchanger. 

The refrigeration system is shown schematically in figure 1 . 
Liquid freon is taken at low pressure and temperature from 
the bottom of the flash cooler and is pumped through a series 
of pipes (approximately 240 ft long) and headers leading 
to tunnel heat exchanger coils. At the design flow rate of 
150 lb/sec, it takes the liquid freon approximately 50 sec to 
travel from the flash cooler to the heat exchanger. The heat 
exchanger, located inside the tunnel, consists of 80 individual 
heat exchanger panels arranged in a pleated design. Rectangular 
sheet fins attached to copper coils carry heat from the tunnel 
air into the copper coils, causing the freon flowing inside the 
coils to partially vaporize before leaving the heat exchanger 
and returning to the flash cooler. 

In the flash cooler the liquid freon drops to the bottom of 
the tank to be again recirculated by the pump. The freon vapor 
in the flash cooler flows through damper valves into an 
accumulator tank. The damper valves control the tunnel 
temperature by regulating the flow in the compressor loop. 
From one to six parallel-operated damper valves can be used 
depending on tunnel heat load. As the freon vapor leaves the 
accumulator, it enters the compressor where its pressure and 
temperature are raised. A super-heated freon vapor then leaves 
the compressors and enters the tube-in-shell condensers where 
cooling-tower water is used to remove heat from the high- 
pressure, high-temperature freon vapor. The heat removed 
from the freon causes the refrigerant vapor to condense into 
a saturated liquid. 
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Figure I .—Refrigeration system for Icing Research Tunnel. (Compressor-condenscr-economizer unit shown inside dashed lines.) 


Freon from the condenser may flow into a two-stage 
economizer (compressor interstage cooler) or may bypass the 
economizer and flow back to the flash cooler depending on 
the operating conditions. When the economizers are used, 
liquid freon enters through a set of orifices which results 
in a pressure drop and causes part of the fluid to flash to the 
vapor state. The flash vapor is routed back to the inlet of the 
second and third stages of the centrifugal compressors. The 
liquid part of the freon leaves the economizers at a low pressure 
and temperature and returns to the flash cooler where part of 
it flashes to vapor and again becomes available to be drawn 
through the damper valves for a restart of the vapor cycle. 
The liquid part of the refrigerant in the flash cooler is 
available to be pumped to the tunnel heat exchanger for tunnel 
heat removal. 

The four-stage centrifugal compressor, a two-pass tube-in- 
shell water-cooled condenser, and a two-stage economizer 
(used to improve the efficiency of the process) are all part of 
a self-contained unit. One compressor-condenser-economizer 
unit is shown schematically within the dashed lines of figure I . 
A 1500-hp induction motor drives each unit. At present the 
refrigeration system contains 13 such units which can be used 
during a run. 


Computer Model 

Model Description 

A computer model was made of the IRT Refrigeration 
System and was used to formulate and test the IRT temperature 
control design presented here. This model can be used for 
further automatic control studies in areas such as startup and 
shutdown. 

The computer model uses algebraic and time differential 
equations to describe steady-state as well as transient con- 
ditions. A lumped -parameter approximation of the system 
components is made using control volumes. Transient equa- 
tions of conservation of mass, conservation of energy, and 
equations of state are used to solve for the temperature and 
pressure for these components. Appropriate flow equations for 
the vapor flow through the damper valves (control valves), 
condensation flow in the condenser, and evaporative flow in 
the flash cooler and economizer are used. Line pressure drops, 
and transients that have no effect on the control design were 
neglected. A schematic of the computer model is shown in 
figure 2. A detailed description of the refrigeration system 
computer model with equations is presented in appendix A. 
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Figure 2 — Simulation diagram for Icing Research Tunnel. (Compressor- 
condenser-economizer unit shown inside dashed lines.) 


The computer simulation was programmed in the CSMP3 
simulation language. This simulation language is well suited 
for modeling systems that can be described by ordinary time- 
dependent nonlinear differential equations. 

Model Validation 

The refrigeration system was instrumented for pressure and 
temperature at various points of interest primarily to obtain 
a valid compressor characteristic curve at the 6650-rpm design 
operating point. These temperature and pressure data were 
also compared with the predicted steady-state results of the 
computer model. 

Table I compares data from a tunnel run with those from 
a similar computer run. The two runs closely agree for all 


TABLE I.— TUNNEL DATA COMPARED WITH COMPUTER RESULTS 


Variable 

Tunnel run data 

Computer results 

Temperatures, °R 



Tunnel 

460 

460.0 

Flash cooler 

448 

454.6 

Compressor inlet 

453 

454.5 

Compressor exit 

679 

676.7 

Water exit 

529 

530.8 

Pressures, psia 



Flash cooler 

19.3 

21.2 

Compressor inlet 

8.3 

8.8 

Compressor exit 

99.0 

98.1 

Compressor motor current, A 

250 

175.2 

Compressor loop flow, lb/sec 

a 24 

23.6 


variables except the compressor motor current. A much higher 
motor current was needed in the actual tunnel run than was 
predicted by the computer model. This indicates that the 
compressor-motor unit was probably operating at a much 
lower efficiency (40 percent) than was used in the simula- 
tion (70 percent). This was not surprising since the system 
was operating at off-design conditions. 

Control Development 

System Characteristics 

The IRT Refrigeration System is not well suited for 
automatic control mainly because of the numerous operating 
constraints of the system: 

(1) Operating constraint caused by compressor surge 

(2) Upper limit of compressor inlet (accumulator) pressure 
of 12 psia (because of seal limitations) 

(3) Upper limit of condenser pressure of 125 psia (because 
of structural limitations) 

(4) Upper limit of compressor motor current of 300 A 

Compressor surge and compressor inlet pressure are the 
most frequently encountered constraints. They can be exceeded 
whenever the tunnel operating temperature is changed too 
rapidly. The compressor can surge when the tunnel temper- 
ature is increased, and the compressor inlet pressure can 
exceed 12 psia when the tunnel temperature is decreased. This 
study focused primarily on the compressor inlet pressure 
constraint because the compressor surge was not well defined 
and will not be as much of a problem. 

Controller Description 

The controller was designed to bring the tunnel air to the 
desired operating temperature as fast as possible without 
exceeding any system constraints. This was best accomplished 
by changing the temperature at a constant rate. Moreover, to 
minimize time delays, the controller was designed to bring 
the temperature to its new steady-state value with little over- 
shoot and minimum settling time. 

Standard control analysis techniques were not used in this 
study because the system is nonlinear. Instead, various temper- 
ature control schemes were tried on the computer. In the scheme 
that was selected, the outer control loop controlled the tunnel 
temperature, an intermediate loop regulated the flash cooler 
temperature, and an inner loop controlled the valve position. 

The flash cooler temperature is used as the intermediate loop 
because it closely approximates the tunnel temperature in 
steady state and responds immediately to a change in control 
valve position. The tunnel temperature responds more slowly 
because the large thermal capacity of the tunnel heat exchanger 
introduces considerable temperature lag between the flash 
cooler and tunnel temperatures. 

In addition to the thermal lag, there is a temperature trans- 
port lag caused by the fluid flowing into the heat exchanger. 


Calculated from data. 
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Figure 3.— Tunnel temperature controller. 


It takes 50 sec for the liquid freon to travel from the flash 
cooler to the heat exchanger. This long dead time, if it were 
a predominant effect, would be a serious hindrance in the 
design of an effective automatic temperature controller. 
However this is not the case here because any immediate freon 
temperature change is due to a change in control valve position. 
The valve position, which affects the flash cooler pressure and 
hence flash cooler temperature, also affects the heat exchanger 
temperature because both pressures are essentially the same. 
Therefore the temperature of the freon in the heat exchanger 
changes at the same time as the flash cooler temperature and 
is not dependent on the fluid flowing in from the flash cooler. 
The temperature of this fluid should have only a minor effect 
on the temperature transients in the heat exchanger. The 
computer model therefore neglects this transport delay. 

The setpoint for the flash cooler temperature control is 
generated by the outer control loop. At steady state, the flash 
cooler temperature generally is 5° lower than the tunnel 
temperature. Therefore to account for this difference, the flash 
cooler temperature setpoint was made 5° lower than the tunnel 
temperature setpoint. Because the difference is not quite constant 
and may change somewhat for different operating condi- 
tions, a proportional temperature control was added. The 
proportional controller not only corrects for the steady -state 
temperature difference, but also speeds up the control response. 

Because the flash cooler fluid is saturated, its vapor pressure 
is strictly a function of temperature. As a result, the flash 
cooler temperature control can be converted to a pressure con- 
trol. This conversion should improve the controller response 


because the flash cooler pressure response is faster than the tem- 
perature response. Also, pressure control may be more desirable 
than temperature control, from a measurement perspective. 

Valve position was added as an inner control loop so that 
the system could be switched to manual control. Since the 
dynamics of the inner loop are much faster than those of the 
other components in the temperature controller, the inner loop 
has essentially no effect on the overall response. 

A block diagram of the Tunnel Temperature Controller is 
presented in figure 3. The AT block represents a gain constant 
for the outer temperature control loop. The KS block repre- 
sents a correction for the steady-state temperature differences 
between the flash cooler temperature and the tunnel temper- 
ature. The “sat curve” block represents a lookup table which 
converts freon temperature to pressure. The “LIMIT” block 
limits the flash cooler pressure error signal before it enters 
the Control Algorithm block. The block output M k represents 
the control valve position setpoint. The KV block is a gain 
constant for the valve position control loop. The equations for 
the control algorithm of the pressure control loop follow. 

Control Algorithm 

The control algorithm is depicted in the block diagram shown 
in figure 3 as follows: 

M k ~2 = M k - j 
M k -x=M k 
er k - 2 = er k - } 

er k - 1 = er k 

Test for setpoint step change; if step, reset 

er k _ i = er k _ 2 = 0 and A /*- 1 = M k + FM 

Then continue 

FM= GO -0.799 er k _ { + 1.998 er> : _ 2 ) 

M k = 1.3 M k _ { - 0.3M*_ 2 + Ger k 4- FM 

where G is a controller gain which is variable and a function 
of tunnel heat load, the control valve position, and the number 
of compressors in operation; er is the control algorithm input; 
M is the control algorithm output; the subscript k-n represents 
the (k-n) xh sampling instant. The algorithm of the pressure 
control loop was not tuned for tunnel pressure response, but 
rather for satisfactory tunnel temperature response. 

Control Run Results 

The controller was tested at various operating conditions 
and for various setpoint changes. When the controller is 
implemented in the field, the tunnel temperature changes will 
probably be made with step changes in the temperature 
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TABLE IT.— NOMINAL OPERATING CONDITIONS 
[Step setpoint change.] 


Initial temperature, °R 490 

Tunnel heat load, Btu/sec 1200 

Compressor used 1 

Control valves used 3 

Evaporation constant 10.0 

Control constants 

Error limit 8.0 

Sample rate, samples/sec 2 

Valve loop gain, KV. 10.0 

Outer loop gain, AT 15.0 


setpoint. Consequently, all but one of the computer runs that 
are presented are for step setpoint changes of various sizes 
and at various operating temperatures. 

The evaporation rate in the flash cooler and the heat transfer 
coefficient in the tunnel heat exchanger were two parameters 
whose values were considered to be very uncertain. Since these 
parameters were deemed important to the system response, 
control runs were made where these parameters were varied. 
Other important system parameters were also investigated: the 
tunnel heat load, the number of compressor units in operation, 
the number of control valve units being used, the effect of 
controller gain (change in the limit of the pressure error signal), 
and controller sampling rate. Because of the uncertainty of 
the compressor characteristic curve, a control run was made 
where the slope of the compressor curve was increased by 100 
percent to determine its effect on the closed-loop response. 

A typical tunnel operating condition and controller setting 
were selected and termed as the nominal case. Nominal 
operating parameters are presented in table II. The time 
response for runs with different operating conditions are 
compared with this nominal case whenever applicable. 

Run 1— Nominal Case 


The first run presented is for the nominal case with the tunnel 
operating at a temperature setpoint of 490 °R. At 80 sec after 



Time, sec 

Figure 4.— Tunnel temperature response for setpoint decrease from 490 
to 450 °R. 


the start of the run, the tunnel temperature setpoint was stepped 
from 490 to 450 °R. The controller performed according to 
design expectations in that the tunnel temperature response 
(fig. 4) was similar to a ramp and (very importantly) the 
response showed no overshoot. It took 485 sec (from 80th to 
565th sec) for the tunnel temperature to reach its new setpoint 
of 450 °R. The control valve movement for the run is pre- 
sented in figure 5. The valve started from an initial position 
of approximately 2° and moved to a final position of approxi- 
mately 4.7°. The compressor inlet (accumulator) pressure is 
a critical parameter in these runs because its upper limit of 



Figure 5.— Control valve movement for tunnel temperature setpoint decrease 
from 490 to 450 °R. 



Figure 6.— Accumulator pressure for tunnel temperature setpoint 
decrease from 490 to 450 °R. 



Figure 7.— Flash cooler pressure change for tunnel temperature setpoint 
decrease from 490 to 450 15 R. 
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12 psia can often be exceeded during the transients. The 
accumulator pressure (fig. 6) was initially 8.1 psia. At the 
start of the transients it reached a peak of 10.6 psia and then 
settled to a new steady-state value of 7. 1 psia. The 12-psia limit 
was not exceeded during this run. The flash cooler pressure 
is provided in figure 7 as a matter of interest. 

Run 2— Compressor Change 

Run 2 was identical to run 1 except that two compressors 
were used instead of one. This change increased the tunnel 
temperature response time by about 45 sec (fig. 8). The 
temperature reached steady state at about 520 sec with a slight 
overshoot. (In run 1, there was no overshoot.) As was expected, 
the system responded faster when more compressors were 
added to the system. 

Run 3 — Control Valve Change 

Run 3 was also identical to run 1 except that four control 
valve units were used instead of three. The tunnel temperature 
(fig. 9) reached steady state at approximately 450 sec with no 
overshoot. The response is at least 100 sec faster than the 
nominal case (fig. 4) of run 1. 


Run 4— Change in Operating Point 

Because the system is nonlinear, run 4 was made to show 
controller performance at other operating points. Run 4 started 



Time, sec 

Figure 8.— Tunnel temperature response for tunnel temperature setpoint 
decrease from 490 to 450 °R with two compressors. 



Time, sec 


Figure 9.— Tunnel temperature response for tunnel temperature setpoint 
decrease from 490 to 450 °R with four damper valves. 



Figure 10.— Tunnel temperature response for tunnel temperature setpoint 
decrease from 520 to 470 °R. 



Figure II .—Accumulator pressure change for tunnel temperature setpoint 
decrease from 520 to 470 °R. 



Figure 12.— Accumulator pressure change for tunnel temperature setpoint 
decrease from 520 to 470 °R with controller gain reduced by 50 percent. 


from an initial operating point of 520 °R. At 80 sec the setpoint 
was stepped from 520 to 470 °R. The tunnel temperature 
response (fig. 10) showed a slight overshoot. Even though the 
setpoint change was 10° larger for this run than for run 1, 
the response is more than 140 sec faster because the required 
change in cooling power is smaller at the higher temperatures. 
One other observation is noteworthy. The higher tunnel 
operating temperature generally raised the pressure and 
temperature of the entire refrigeration system, increasing the 
likelihood of the accumulator pressure limit being exceeded 
during transients. Figure 1 1 shows the accumulator pressure. 
It exceeds its 12 psia limit and briefly rises above 17 psia. 
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Figure 13.—' Tunnel temperature response for tunnel temperature setpoint 
decrease from 520 to 470 °R with controller gain reduced by 50 percent. 



Figure 14.— Tunnel temperature response for tunnel temperature setpoint 
ramped from 520 to 470 °R in 100 sec. 



Figure 15. — Tunnel temperature response for tunnel temperature setpoint 
decrease from 460 to 440 a R. 


this facility is automated. The new accumulator pressure limit 
will then be approximately 20 psia. 

Run 6— Setpoint Ramp 

In the previous runs all temperature setpoint changes were 
step changes. To test the controller for a ramp setpoint change, 
run 6 was made with the tunnel temperature setpoint ramped 
from 520 to 480 °R in 100 sec. The tunnel temperature response 
(fig. 14) is similar to a step setpoint change. 

Run 7— Small Setpoint Change (Actual Tunnel Run 
Comparison) 

A small setpoint change was made from 460 to 440 R at 
80 sec. It takes 380 sec or 6.3 min for the temperature to reach 
steady state at its new operating point (fig. 15). This run was 
compared with data from an actual run where the operator 
controlled the tunnel temperature manually. The system 
operator took 20 min to manually change the tunnel tempera- 
ture from 460 to 440 °R. This is more than three times longer 
than by automatic control. 

Run 8— Heat Load Ramp 

The tunnel fan speed is often changed during tunnel operation 
to accommodate wind tunnel test requirements. Because 
the fan is the major tunnel heat source, changing the fan speed 
also changes the tunnel heat load. The effect of this heat 
load change on the temperature control was studied by ramping 
the heat load from 800 Btu/sec (at 100 sec) to 1500 Btu/sec 
(at 200 sec). The setpoint was held at 460 °R during the ramp. 
The automatic controller was able to maintain the tunnel 
temperature near its setpoint with essentially no temperature 
change (fig. 16). Figure 17 shows the valve position change 


? 460.2 r- 
E 

| 460.0 
8 . 

£ 459.8 


800 Btu/sec 



1500 Btu/sec 


100 200 300 400 500 

Time, sec 


600 


Figure 16. — Tunnel temperature variation during a tunnel heat load increase 
of 700 Btu/sec. 


Run 5— Controller Gain Change 

A lower controller gain is needed to keep the accumulator 
pressure in run 4 below its limit. Consequently, run 4 was 
repeated with the controller gain (limit on the pressure error 
signal) lowered by 50 percent. The reduced gain kept the 
accumulator pressure transient below the 12-psia limit 
(fig. 12), but it also produced a slower temperature response 
(fig. 13; about 340 sec slower than in run 4). The gain will 
probably not have to be lowered, for this case, because a 
refurbished set of compressors should be installed by the time 



Figure 17.— Control valve position during a tunnel heat load increase of 
700 Btu/sec. 
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made by the controller in order to maintain the temperature 
at its setpoint during the heat load ramp. The valve moved 
from an initial position of approximately 2.4° to a final posi- 
tion of 4.9°. 

Run 9— Setpoint Increase 

The general operating procedure for all tunnel runs is to 
start at ambient conditions and then decrease the tunnel 
temperature to the desired operating point. Therefore, all of 
the computer runs presented thus far were for decreases in 
the temperature setpoint. Occasionally, however, the tunnel 
temperature needs to be increased to a new operating point. 
Figure 18 shows such a computer xun where the temperature 
setpoint is increased from 450 to 490 °R at 700 sec into the 
run (run 9). The quality of the response is similar to the step 
decreases shown in the previous runs. For this run it takes 
690 sec for the temperature to go from 450 to 490 °R. There 
is no overshoot. During this excursion the accumulator pres- 
sure transient moves away from the 12-psia limit. This is 
typical for all runs where the tunnel temperature is increased. 
Because the accumulator pressure is not a problem for these 
runs, it is not examined when the tunnel temperature is 
increased. Compressor surge, however, can be a problem for 
large, rapid tunnel temperature increases. The problem was 
investigated; but because the compressor surge line for the 
present compressors is not well defined, the run results from 
this investigation are not presented here. If compressor surge 
becomes a problem for the automated system, the controller 
gain will have to be reduced for certain runs involving tunnel 
temperature increases. 

The response to a setpoint increase in the tunnel temperature 
is often limited by the present configuration of the refrigeration 
system. For this system, the control valve usually needs to 
be operated near its closed position— a practice that restricts 
any valve movement in the closed direction and, consequently, 
limits the response for tunnel temperature increases. The valve 
movement for run 9 is shown in figure 19. Prior to the setpoint 



Time, sec 

Figure 18.— Tunnel temperature response for tunnel temperature setpoint 
increase from 450 to 490 °R. 



Time, sec 

Figure 19.— Control valve position for tunnel temperature setpoint increase 
from 450 to 490 °R. 

increase, the valve position was 4.6°, As soon as the setpoint 
increased, the valve nearly closed and remained so until the 
tunnel temperature had almost reached the new setpoint value. 
At that time, the valve moved to a final steady-state position 
of 1.4°. 

Runs 10 to 13 — Additional Runs With Setpoint Increases 

Other runs tested different temperature setpoint increases 
similar to that in run 9. The tunnel temperature response for 
each of these runs was very satisfactory with essentially no 
overshoot or oscillation as is shown in figure 20. The operating 
conditions for these runs were similar to the previously 
presented runs. 

Small Setpoint Disturbances 

The next set of runs were made to further demonstrate the 
robustness of the controller and to show the effect of various 
parameters on closed-loop performance. For these runs the 
setpoint was stepped from 460 to 440 °R. Run 7 will be 
considered as the reference or nominal case, and all parametric 
variations of the controller and the system will be compared 
to figure 15 of run 7. 

Controller Variation Runs 

Figure 21 shows the effect of two important controller 
settings on the temperature response. First, the controller 
sampling rate was reduced from 2 samples per sec to 1 sample 
per sec (the limit on the error signal was doubled to offset 
the gain change caused by the sampling rate). The sampling 
rate change showed virtually no effect on the system. Next, 
the limiter output of the error signal was changed from its 
nominal value of 8 to a value of 12 (see fig. 3). The tem- 
perature shows no overshoot and reaches its setpoint value 
about 80 sec sooner and its steady-state value almost 150 sec 
sooner than for the nominal gain setting. 
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(a) Setpoint increase from 440 to 480 °R at nominal operating conditions. (c) Setpoint increase from 440 to 460 R with three compressors and a tunnel 

(b) Setpoint increase from 440 to 460 °R with controller gain at 150 percent. heat load of 2000 Btu/sec. 

(d) Setpoint increase from 450 to 490 °R with tunnel heat load at 2000 Btu/sec. 

Figure 20.— Tunnel temperature response for various tunnel temperature setpoint increases. 


System Variation Runs 

Two system parameters which were thought to be important, 
and yet whose values are very uncertain, are the evaporation 
rate in the flash cooler and the heat transfer coefficient in the 
tunnel heat exchanger. A large part of the heat transfer 
coefficient uncertainty comes from the approximate calculation 
of the fluid quality in the heat exchanger. Figure 22 shows 
the temperature response for a nominal value of the 
evaporation rate constant, KVE, and for a 100 percent increase 
of KVE. The increase had little effect on the temperature 
response. The temperature response increased slightly when 



Figure 21.— Tunnel temperature response for various controller settings 
(sampling rate and gain). 


the heat transfer coefficient was changed to twice its nominal 
value, but the results are not presented here. 


Heat Load Variation Runs 

The tunnel heat load is mostly a function of tunnel fan speed 
(rpm) and, to some extent, the outside ambient temperature. 
Figure 23 shows the controller performance for different tunnel 
heat loads— 1200 (nominal), 2000, and 800 Btu/sec. As would 
be expected, the lowest heat load shows the best response. 



Figure 22.— Tunnel temperature response to variations in evaporation rate. 
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Tunnel 



Figure 23.— Tunnel temperature response for various tunnel heat loads. 


Compressor Variation Runs 

Because of the uncertainty in the compressor characteristic 
curve, and the likelihood that a refurbished compressor will 
be installed in the system, a control run was made where the 
slope of the current compressor characteristic curve was 
increased by 100 percent (values for compressor head were 
doubled). The temperature response with the new compressor 
slope is shown in figure 24. The response shows no overshoot. 
The temperature reaches its steady state approximately 30 sec 
sooner than for the run with the nominal compressor charac- 
teristic curve. 

Noise Injection 

The controller was also tested with Gaussian noise added 
to the control signals of the flash cooler and to the tunnel 
temperature. The added noise was considerably higher than 
can ever be expected in the actual instrumentation. Since its 
effect on the controller was essentially negligible, the results 
are not presented here. 


Slope of 
compressor 



Figure 24.— Tunnel temperature response to a change in compressor 
characteristic curve. 


Concluding Remarks 

Although the Icing Research Tunnel (IRT) was designed to 
operate at temperatures below 440 °R, most tunnel runs are 
currently made above 440 °R. Consequently, operating con- 
straints are frequently encountered. These constraints are 

(1) Operating constraint caused by compressor surge 

(2) Upper limit of compressor inlet (accumulator) pressure 
of 12 psia 

(3) Upper limit of condenser pressure of 125 psia 

(4) Upper limit of compressor motor current of 300 A 

The most frequently encountered operating constraint 
appears to be the upper limit of the compressor inlet pressure. 
To avoid this constraint, an additional compressor unit is 
brought on line. This allows the control valves to operate at 
a more closed position, thereby causing a larger pressure drop 
and, consequently, a lower compressor inlet pressure. Since 
this practice brings more compressors on line than would 
otherwise be needed, it reduces the efficiency of the operation, 
and frequently steam must be injected into the tunnel to 
maintain the desired tunnel temperature. The compressor inlet 
pressure also poses a problem for the automatic temperature 
control. The controller must be designed so that this pressure 
constraint is not exceeded during control transients. 

The primary incentive for adding an automatic temperature 
control is to make more efficient use of the tunnel run time. 
Consequently, the automatic temperature control was designed 
to reach the desired operating point as quickly as possible 
without creating excessive transients that would exceed system 
constraints. The objective was to ramp the tunnel temperature 
to its desired value with a constant rate and little overshoot. 
For the runs that were made in this study, the control design 
met this objective. Study results showed that the automatic 
temperature controller can change the tunnel temperature in 
less than half the time needed by manual control. We compared 
a computer run with an actual, manually controlled tunnel run. 
For the computer run, the automatic control lowered the tunnel 
temperature from 460 to 440 °R in only 6.3 min, whereas 
for the actual run with manual control, 20 min were needed 
to make the same change. The automatic controller, however, 
frequently caused a transient, excessive compressor inlet pres- 
sure when the controller gain was set too high. 

System constraints could be avoided if automatic control 
were used for selected runs only— when there is no danger 
of exceeding system operating limits. The alternative would 
be to use automatic control at all times but with a low controller 
gain. This would avoid the large transients which exceed the 
operating limits. However, with a low gain, the system would 
respond slowly, and the tunnel run time would not be 
significantly reduced compared with manual operation. The 
only benefit of using an automatic control with the low gain 
option, would be to free the system operator from the task 
of manually controlling the tunnel temperature. 
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In the interest of saving tunnel run time and making the 
compressor operation more efficient, automatic temperature 
control should be added to the facility, and the facility should 
be modified to allow a higher compressor inlet pressure. 
Currently the compressors are being overhauled to make them 
more suitable for higher tunnel temperature runs. The over- 
hauled compressors will permit operation at a suction header 
(accumulator) pressure of up to 20 psia. No major change in 
the control should be required for operation with the 
overhauled compressors. 

The automatic control study presented here does not address 
an automatic control design for startup and shutdown or any 
decision-making capabilities whereby the computer would take 
the necessary action to avoid the various operating constraints. 


These decisions and tasks will still have to be performed by 
the system operator. It is recommended, however, that, after 
successfully implementing this initial phase of the automatic 
control, automation of these remaining tasks should also be 
considered. The present computer model was made with 
sufficient detail to allow these additional control studies to be 
performed with minimal changes. 


Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, August 13,1990 
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Cooling tower water 


Appendix A 

Computer Simulation Equations of the IRT Refrigeration System 

Refrigeration Cycle 

The refrigeration cycle for the compressor loop is depicted in 
figure 25. Saturated gas in the flash cooler (evaporator) moves 
through the control valves to the accumulator, or compressor 
inlet. The gas is then compressed to the superheated region 
at the compressor outlet. From the compressor outlet the gas 
flows into the condenser region, where heat is removed and 
the gas condenses to a saturated liquid. The liquid then 
flows into the economizer region (fig. 25 depicts a two-stage 
economizer), where it loses more heat at a reduced pressure 
and temperature. From there it flows back to the flash cooler. 

A diagram of the computer model is shown in figure 26. 

Basic Computer Model Equations 

The following basic fluid flow and energy equations are used 
repeatedly throughout the simulation. The equations are based 
on a control volume concept of finite difference equations, 
where i and o represent flow into and out of the control 
volume, respectively. 

Continuity equation 


dM 

Hi 


= W; - W n 


( 1 ) 


Energy equation 


dE 


= W t H t - Wjr o + Q 
dt 


( 2 ) 



and Q represents external heat conducted into or out of the 
control volume. For convective heat transfer, 


Q — hS(Tft m6 — T wan ) 


(3) 


where the total enthalpy E is defined as the mass M times the 
enthalpy H , or 

E = MH 



where 

h heat transfer coefficient 
S heat transfer area 

Since dH = c p dT , and if the term dM/dt can be assumed to 
be negligible, an alternate form of equation (2) is 


dT 

~dt 


- WT + 


M 


Perfect gas law 


P = pRT 


Orifice equation 


W = 0.6A 



2g c P i (P i ~ P 0 ) 

ZRT { 


(2a) 


(4) 


(5) 
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Evaporation equation 1 

W = KVE ( P sal - P G ) (6) 

where 

KVE = KrA\f — 

V 2itRT L 

The preceding equations are applied in the mathematical 
description of the system components presented next. 


Likewise, the freon gas temperature can be obtained from the 
saturation curve of the PH diagram 

T G , e =f(H Gie ) 

Also, the enthalpy of the evaporating freon can be obtained 
from the freon saturation curve: 

Hv,e =f(T L e) 

And again, since the pressure of the freon in the heat exchanger 
and the flash cooler are essentially the same, 


Flash Cooler Equations (Evaporator) 

Energy equation.— From the energy equation (2), we have 
the following differential equation for the flash cooler 
(evaporator) liquid, where N c accounts for the number of 
compressor units in operation: 


H Gx — H Ve 

Continuity equation. — From the continuity equation (1), we 
have for the liquid 


d(M Le ) 

dt 


= w Ltec + w l>x - w ar - w p 


= N c W Lec H Uec + W u H LiX 

dt 

- W p H u - W C l H Vx 

Since W p = W Lx + W Gx and H Lx = H Le because the 
pressure of the freon in the heat exchanger and the flash cooler 
are essentially the same, we can simplify the above equation 
as follows: 


(i ^ L = N c W Le< H Lec - W c , e H v , e ~ W Cx H U( 

dt 


H Le = 


-*L,e 


Mi 


L,e 


Likewise, the energy equation for the flash cooler gas is 


and for the gas we have 


d(M Gie ) 

dt 


+ W Gre - W M 


Equation of state. — From the definition of density we have 
for the density of the gas p G>e = M G e fV G e , and from the 
density of the liquid, V Le — M Le /p Le . The flash cooler 
consists of both gaseous and liquid freon. Therefore the total 
flash cooler volume V e = V Ge + V L>e . Using these three 
relations, we can now solve for the density of the flash cooler 
gas as follows: 


PG,e — 


M G , e 


M L> e 

PL,e 


d{E Ge ) 

dt 


= Wr.rH 


G,x n G,x 


KcH C ,e + W G .eHv,e 


From the perfect gas law we can solve for the freon gas 
pressure 


From the saturation curve on the pressure-enthalpy (PH) 
diagram for freon (fig. 25), the freon liquid temperature is 
given by 

T Lr e=AH L J 


P G ,e = ZpG,eRT G>e 

The density of the freon liquid varies as a function of the 
temperature as follows: 

p Le = 138.84 - 0.1057^ 

Evaporation rate. — Equation (6) is used to solve for the 
evaporation rate of freon in the flash cooler: 


’See Derivation of Evaporation Equation, p. 19. 


W G , e = KVE e (P^ e -P G ,e) 
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where 


KVE e = KrA e 


V 2wRT Le 


Let us now solve for the mass of the two-phase fluid inside 
the heat exchanger. Assume that the mass of the gas in the 
heat exchanger increases linearly as the fluid flows through 
the heat exchanger. We can then express the differential mass 
along the heat exchanger as follows: 


The freon saturation pressure P sal e can be obtained from the 
saturation curve 

= /( T Le ) 

Heat Exchanger Equations 

Energy equation for freon.— The flow throughout the heat 
exchanger loop is assumed to be equal to the pump flow in 
the loop. The energy equation for the enthalpy at the heat 
exchanger exit can then be written as follows: 

d(E x ) 

= W p (H L>e - H x ) + Q x 



where M x is the mass of the fluid inside the heat exchanger 
and W p is the pump flow rate. The heat flow from the heat 
exchanger wall into the freon liquid, Q x , is shown in the 
following equation: 

Qx = h x S x (T m - T l x ) 

Since freon in the heat exchanger is assumed to be at the 
saturation pressure of the flash cooler, we have 

Ti iX = T L e p Gx = p G e H C x = H V e 

In a two-phase system, X is called the quality of the mixture 
and is defined as the weight fraction of the material in the vapor 
phase, or 



M 


From the definition of X, we can likewise define the enthalpy 
of the gas and liquid mixture H x as 

H x = (l — X)H Lx + XH Gx 


d{M G ) — Pg a g d? 
d{M L )= PL {A x ~A G )dt 

where 

A g gas flow area 
A x total flow area 
dt differential length 

Since it can be assumed that the gas density is constant along 
the heat exchanger, we can state that the cross-sectional flow 
area of the gas for the heat exchanger of length L varies linearly 
from zero at the entrance, to some value A Co at the exit. 
Then at any length f of the heat exchanger, 



Now the total mass of the gas in the heat exchanger is 
M c =J^(dM c ) = f p c A G d( 

J 0 


= PgAc A l 

L Jo 


Similarly the mass of the liquid in the heat exchanger is 


A h = P L l(a x - 

The total mass of the fluid in the heat exchanger is 


Rearranging the above equation we can solve for the quality 
of the fluid in the heat exchanger as follows: 


M x — + M g — PiA x L + A Go 



(. Pg — Pl) 


where 


X r = 




AH X = H a>x - H L>X 


From the definition of quality, the quality at the heat exchanger 
outlet is 


X fi = 


dM { 


G.o 


dM Go + dM Lo 


Pg a g,q 

PG A G,o + Pl( A x ~ A G,o ) 
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Solving for A G , 0 we have 


where 


A G ,o 


XqPjAx 

P C ( 1 - X 0 ) + Xj) L 


After substituting for A G>0 in the previous equation, the mass 
of freon in the heat exchanger as a function of quality is 



( X o\ 

x. 

plK 

K v 

T Pt 


p c ( 1 - X 0 ) + X 0 p L 


where 

V x = AJL 

The following heat transfer coefficient was used for the freon 
side of the heat exchanger (ref. 1): 


K 


0.0225 



W„D X 

AxVx 


0.75 


— (/) 
8 


xjAH x \ 
2 \L X ) 


0.375 


Heat transfer to the fluid is assumed negligible as it flows from 
the flash cooler to the heat exchanger. Therefore, 


T _ + T aJ 

a 2 

the following heat transfer coefficient (ref. 2) was used for 
the air side of the heat exchanger: 


K 


0.1376fc g fWJh \° 61 

D a \ ^-aPa / 


Energy equation for the tunnel air. The tunnel air is 
divided into two control volumes. One volume V a x consists 
of the heat exchanger portion of the tunnel air. The other 
volume V aJ consists of the remaining tunnel air which 
includes the tunnel fan where the heat is generated; T a j is the 
air temperature at the heat exchanger entrance, and T ax is the 
air temperature at the heat exchanger exit. The flow rate in 
the tunnel circuit, W a , is assumed to be the same in both 
control volumes. The energy equation for the heat exchanger 
control volume can be written as follows: 


d(T a>x ) 

dt 


W a (T aJ -T atX )-f- 

Cp.a 


M„ 


H L , x = H L , e 

A H x = A H e 

from the definition of quality we can solve for the gas flow rate: 


where 


M ax Pa^a,x 

likewise the energy equation for the fan portion of the control 
volume can be written as follows: 


W GiX = X x W p 

Since the pump flow is the total heat exchanger flow, we have 
for the liquid flow 


W a {T a<x - T a j) + Q- 
d(T aJ ) C J± 

dt M aJ 


W LtX = W p - W G>X 

Energy equation for the heat exchanger wall. Consider 
the energy balance for the heat exchanger wall. Heat Q a 
flows from the tunnel air into the wall and heat Q x flows out 
of the wall into the freon heat exchanger fluid. The energy 
equation (2a), has the following form for the wall: 

d(T m ) = (Qa - Qx) 

dt MfnCpm 

the equation for the heat flowing from the tunnel air into the 
wall can be written as follows: 

Qa ~ ha^a (T a ~ ) 


where 

M a J ~ PaVaJ 

The heat rate generated by the fan Qf depends on the fan 
speed. 

Pump Flow Loop Equations 

Initially this part of the simulation was more detailed but 
was later simplified because of some occasional computational 
problems and because it was found that a small variation in 
pump flow had little effect on the system. The pump flow was 
therefore made constant, where 

W p = 150 
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Accumulator Equations 

Continuity equation and equation of state.— From the 
continuity equation we have 


d(P M -) _ War - W c 
dt V ac 


Tac 


Compressor Equations 


The pressure in the accumulator can now be obtained from 
the perfect gas law 


^ac PacZRT a , 


Orifice equation.— The orifice equation is used to solve for 
the flow through each control valve unit, or 


fVi 


= Q.6A s f\J ■ 


1j£gAPg± - Pa, ) 

ZRT„, 


where T acJ the temperature in the accumulator, is assumed 
to be equal to the temperature of the freon gas in the flash 
cooler. The effective control valve flow area A v is shown in 
figure 27. The flow into the accumulator is 


Kc = A JyWy 

w'here N v is the number of damper valve units in operation. 


Compressor characteristic curve. —The compressor charac- 
teristic curve shown in figure 28 was obtained from a limited 
amount of available tunnel run data. The ordinate represents 
the compressor head in thousands of feet and the abscissa is 
the compressor volumetric flow rate in thousands of cubic feet 
per minute. 


Compressor flow equations.— The equations are 


HEAD = ZRT^ 




- 1 


Q = /(HEAD) 
ii/ __ QPac 

rr ci 

60 

The compressor outlet flow is equal to the compressor inlet 
flow plus the vapor flow from the economizer, or 


Energy equations . The energy equation for the accumu- 
lator is 

d(E m ) 

-^=W«-H G .e-N c W d H„) 

where 

W Ui . flow into the accumulator 
N c W d total flow leaving the accumulator 


Ko = W ci + W G>ec 

Line pressure drops from the accumulator to the compressor 
inlet and from the compressor outlet to the condenser are 
neglected, or 




Figure 27. — Control valve area. 


There is also no heat transfer from the accumulator up to the 
compressor inlet, or 

T = T 

1 a 1 at 


18x10 3 



Figure 28.— Compressor characteristic curve; speed, 6650 rpm. 
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Compressor energy equations .— The compressor is divided 
into two segments. All of the economizer gas is assumed to 
discharge into the compressor at the entrance of segment 2. 
The equation below is the temperature rise of the fluid for 
compressor segment 1: 


7-1 



where the flow into the condenser is assumed to be the flow 
out of the compressor W co . Also the enthalpy H co out of the 
compressor is assumed to be the same as the enthalpy flowing 
into the condenser. The variables W con and Hi c on are, 
respectively, the flow and the enthalpy of the liquid freon 
leaving the condenser. 

By neglecting the heat capacity of the condenser piping, the 
heat rate removed from the condenser by the water, Q c on > can 
be written as follows: 


Since <W actua i = Vv a {dH iscnxwp[c ) and since dH = c p dT, we 
have for the enthalpy 

H cl =H ac +(^y T cl -T ci ) 

The vapor from the economizer with enthalpy H Vec is 
returned to the compressor and mixes with the compressor 
fluid of segment 1 . From the energy equation, we can express 
the enthalpy of the mixture as follows: 

„ W a H c i + W Gtec H Vi€C 


or the temperature of the mixture 

W a T c] + W Gi€C T L>ec 


Fluid with enthalpy H c2 and temperature T cl enters the second 
segment of the compressor and rises to a temperature T c0 at 
the compressor exit, or 


Ccon — ^con^con(^L,con ^vr) 

The average condenser water temperature T w is calculated as 
follows: 

T w = 0.5(7;, 4- T wo ) 

The heat transfer coefficient for the condenser is too complex 
to be calculated analytically. An overall heat transfer coefficient 
h con o was calculated based on results of known system 
operating conditions. The 0 assigned to variables represents 
tunnel run data which was obtained either directly or calcu- 
lated from data. This data was used to calculate h con Q. 
The equations used in the calculation h con 0 are based on 
the steady-state part of the energy equation (2) and also 
equation (3), or 

Ccon.O = ^con,0^con(^i.,con,0 — ^h>,o) 

But in steady state the heat rejected at the condenser is equal 
to the heat added by the compressor plus the heat added in 
the tunnel, or 

Q CO n,0 = 4.0 + G/.O = W co n fi( H co.O “ H G.e. o) + 4.0 
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The enthalpy at the compressor exit then is 


solving for /i CO n,o we have 


^con.O 


Qjq + (Hco,0 ~ ^G,e,o)^con,0 
^con ( ^L,con ( 0 


H C0 = H c2 + (^j(T CB -T c2 ) 

Condenser Equations 

Energy equations.— The energy equation for the condenser 
can be written as follows: 


A heat transfer coefficient was used which is a function of 
the cube root of condenser flow, or 

l — if w 0.333 

"cort.O — " con,0 

Using the above equation to solve for K h gives 


d(Ei' con) _ 

dt 


= W'con(W G ,con “ " Qo 


K„ = 


“con,0 

iif 0.333 
W con.O 


77^ )C on 


^I,con 


Mj 


L, con 


The value of K h calculated above from tunnel run data is now 
used to calculate h con for all other operating conditions, or 
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h con - KhW° c on 33 . Or by eliminating K h we can write the heat 
transfer coefficient in the condenser as follows: 


Equation of state.— From the definition of density we can 
obtain the density of the gas as follows: 


^con ™ 


con,0 


w 

rr con 


w f 


con.O 


0.333 


PC, con 


Mqq on 
r G,con 


The freon leaving the condenser is assumed to be saturated 
liquid. We can therefore obtain the corresponding temperature 
of the liquid from the saturation curve for freon, or 

^I r con = /(^L,con) 

The energy equation for the water leaving the condenser is 


From the perfect gas law, we can now solve for the pressure 
of the freon inside the condenser, or 

^con = PG, con ZRTc xon 

^sat.con the saturation pressure of the freon in the condenser. 
It can be obtained from the freon saturation curve of the PH 
diagram: 


d(T Mt ) 

dt 


(T wi - T wo ) 


Kpw 



where W w is the water flow rate through the condenser and 
T wo is the outlet water temperature for the condenser. The 
inlet water temperature T wi is assumed to be constant. 

Continuity equation.— From the continuity equation 


<*(M C on) 

dt 


= W„ 


w rt 


^sat.con “ /( con) 

Economizer equations 

The system uses two economizer stages where pan of the freon 
liquid is flashed to gas and returned to the inlet of compressor 
stages 2 and 3. For simplicity, these will be approximated in 
the simulation by one economizer only. Because the fluid 
volumes in the economizers are relatively small, the transients 
in the economizer are neglected. 

Flow equation.— The economizer gas is returned to the 
center stage of the compressor. The average pressure at the 
center stage is estimated as follows, where K e j is a constant: 


where M con is the total mass of freon inside the condenser and 
^con is the liquid freon leaving the condenser (assumed equal 
to the freon condensation rate). 

The relative volume occupied by the gas and liquid inside 
the condenser is assumed to remain constant (i.e,, V Lcon is 
constant). The mass of the gas and the liquid freon inside the 
condenser can then be obtained from the following two 
equations: 

Ml'Con ~ ^L.con Pl 
^G, con ^con ^L,c on 

Freon condensation rate. —Equation (6) is used to solve for 
the condensation rate of freon in the condenser, or 


P e c,j = K ed P C0 + (1 - K ed )P ac 

The frictional loss of the economizer gas discharge flow is 
assumed to be negligible. Consequently, the pressure P ec 
inside the economizer is assumed to be equal to the discharge 
pressure at the compressor P ec j, or 

P — P 

ec r ec,d 

Since the fluid in the economizer is at saturation, it is assumed 
that the economizer saturation pressure P M et . is equal to P ec : 


Since economizer transients are neglected, the steady-state 
continuity equation results in 


W'con A’K£’ t0 n(/’ CHn / > sa , con) 


w L ,ec = *n,,„ - w Ge 


where 


KVE con 



8c 

27rRT Lxou 


Energy equation.— From the steady-state energy equation, 
we can write for the economizer 

^con^Z.con — ^G.ec^G.ec ~ ^L.ec^L.ec = 0 
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Substituting the previous expression for W Lei and after 
rearranging, the freon vapor returning to the compressor is 


W= Ap 



Wa„ = 


^on^.con - H l „.) 
H G ,cc - H L ,ec 


The fluid inside the economizer is assumed to be in the 
saturated region. The freon temperature and the freon enthal- 
pies of the liquid and gas can then be obtained from the 
saturation curves, or 


= /(*W-> 


From gas kinetic theory, the kinetic energy of a gas molecule is 



where T is the temperature and k B is the Boltzmann constant. 
From the perfect gas law PV = nk B T; kinetic energy KE is 
then equal to 


Hu>c=f(T L ,,) 

H GriT =f(T Liec ) 

Derivation of Evaporation Equation (from ref. 3) 

A nonrigorous derivation of the equation used in calculating 
the evaporation flow is given here. Because of the complexity 
of the molecular motion, this equation cannot give good 
quantitative results. The equation is used mainly to account 
for the transient dynamics of evaporation in the flash cooler 
and the condenser. According to gas kinetic theory the number 
of molecules striking a unit area per unit time is 




and from the definition of kinetic energy we have 

Then equating the two kinetic energy equations gives 
(3/2 )P(V/n) = mC 2 / (2g c ) , and substituting for n(m/V) = 
M/V — p = P/{RT) and rearranging we have 


C = yJ~2gJlT 


where 

n/V number of molecules per volume 
C a mean velocity of the molecules 

The mass of the molecules striking a surface area A per unit 
time can be expressed as follows: 



where 

Na Avogadro’s number 
m mass of molecule = mw/Na 
mw molecular weight 

By expressing the gas density p as follows, 


M nm 



we can now rewrite the mass flow' rate W as 


The mean velocity C u is related to the root mean square 
velocity C as follows: 



Now' since W = Ap(C a /4) — A(P/R/T) (C u /4) and after 
substituting for C (l and rearranging we have for the mass 
flow rate 


W = P/l-\ / 

V 2ttRT 
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The preceding equation describes the ideal gas condensation 
rate. It assumes that all of the molecules that strike the surface 
of the liquid will condense. In actuality some of the molecules 
will be reflected back. To account for this, the equation is 
modified as follows: 

W= KrAP\ / 

V 2-kRT 

where the reflection coefficient Kr is the fraction of the 
molecules that hit the surface and are absorbed by the liquid. 
Now suppose we have a closed vessel where the liquid resides. 
Then at equilibrium the condensation rate W con is equal to 
the evaporation rate W ev . Let us suppose further that the 
evaporation flow rate can be expressed similarly except that 
T is the temperature of the liquid and P is the saturation 
pressure of the liquid. We can now express the net flow 


between the gas and the liquid as follows: 



Since the temperature of the gas T G can be assumed to be 
equal to the temperature of the liquid T Lt we then have 
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Appendix B 
Symbols 


A 

flow area, ft 2 

vP c 

freon gas pressure, psf 

A v 

total valve area, ft 2 

P sat 

saturation pressure, psf 

A x 

freon flow area in heat exchanger 

P 

pressure, psf 

C 

root mean square velocity, ft/ sec 

Q 

volumetric flow rate, ft 3 /min 

C a 

mean velocity, ft/sec 

Q 

heat flow rate, Btu/sec 

C P 

specific heat, Btu/lbm- °R 

R 

freon gas constant, 12.776 ft-lbf/lbm- °R 

D 

diameter, ft 

S 

surface area, ft 2 

di 

differential length 

T 

temperature, °R 

€ 

error 

t 

time, sec 

er 

error (limiter output) 

V 

volume, ft 3 

ER con 

fraction of condenser volume occupied by gas 

Kj 

tunnel air volume outside of heat exchanger, ft 3 

/(HEAD) 

compressor characteristic curve 

Va, 

total tunnel air volume, ft 3 

G 

controller gain 

V a , x 

tunnel air volume assigned to heat exchanger, ft 

g 

gravitational constant, ft/sec 2 

K 

freon volume in heat exchanger, ft 3 

gc 

conversion constant, 32.2 lbm-ft/lbf-sec 2 

w 

flow rate, lbm/sec 

H 

enthalpy, Btu/lbm 

X 

quality, dimensionless 

H v 

enthalpy of saturated vapor, Btu/lbm 

z 

gas compressibility factor, dimensionless 

HEAD 

ft-lbf/lbm 

A H 

heat of vaporization, Btu/lbm 

h 

heat transfer coefficient, Btu/sec-fr-°R 

7 

specific heat ratio, dimensionless 

J 

conversion constant, 778 ft-lbf/Btu 

Va 

adiabatic efficiency 

K € M 

flash cooler flow discharge constant 

Vc 

compressor efficiency 

K f 

friction constant 

M 

viscosity, lbm/ft-sec 

KE 

kinetic energy 

P 

density, lbm/ft 3 

Kr 

fraction of molecules absorbed by liquid 



KS 

controller constant 

Subscripts: 

KT 

controller gain constant 

a 

air; adiabatic 

KV 

controller gain constant 

ac 

accumulator 

k 

thermal conductivity, Btu/sec-ft-°R 

c 

compressor 

h 

Boltzmann constant 

ci 

compressor inlet 

L 

length, ft 

CO 

compressor outlet 

t 

variable length of heat exchanger, ft 

con 

condenser 

M 

mass, Ibm 

cl 

compressor segment 1 

M k 

controller output 

c2 

compressor segment 2 

m 

mass of molecule, mw/Na 

d 

discharge 

mw 

molecular weight 

€ 

flash cooler (evaporator) 

N c 

number of compressor units in operation 

ec 

economizer 


pump speed, rpm 

ev 

evaporator 

K 

damper valve units 

f 

fan 

Na 

Avogadro number 

G 

gas 

n 

number of molecules 

i 

inlet 

p 

1 ec,d 

economizer discharge pressure, psf 

k-n 

(£-/i) th sampling constant 


L 

liquid 

w 

water 

m 

heat exchanger metal 

wi 

water in 

o 

outlet 

wo 

water out 

P 

pump 

X 

heat exchanger 

sat 

saturation 

0 

tunnel run data or calculated from tunnel 

V 

valve 


run data 

V 

vapor 
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